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ABSTRACT 

We present new high spectral resolution VLT/UVES spectroscopy and archival 
Cn ! HST/STIS imaging and spectroscopy of the giant Hn region Tol89 in NGC 5398. 

From optical and UV HST images, we find that the star-forming complex as a whole 
contains at least seven young compact massive clusters. We resolve the two brightest 
optical knots, A and B, into five individu al young massiv e clusters along our slit, Al-4 
and Bl respectively. From Starburst99 ijLeitherer et al.f l UV spectral modelling, and 
nebular H/3 equivalent widths in the optical, we derive ages that are consistent with 
the formation of two separate burst events, of ~ 4 ± 1 Myr and < 3 Myr for knots A 
■ (Al-4) and B (Bl), respectively. An LMC metallicity is measured for both knots from 

\^ ' a nebular line analysis, while nebular Hell 4686 is observed in knot B and perhaps 

, in knot A. We detect underlying broad wings on the strongest nebular emission lines 

indicating velocities up to 600 km s . From UV and optical spectroscopy, we esti- 
mate that there are ~95 early WN stars and ~35 ea rly WC stars in Tol89 -A, using 
II i empirical template spectra of LMC WR stars from ICrowther & HadfieldL with the 

Q ■ WC population confined to cluster A2. Remarkably, we also detect a small number of 

| approximately three mid WNs in the smallest (mass) cluster in Tol89-A, A4, whose 

spectral energy output in the UV is entirely dominated by the WN stars. From the 
strength of nebular B./3, we obtain N(O) ~690 and 2800 for knots A and B, respec- 
tively, which implies N(WR)/N(O)~0.2 for knot A. We also employ a complementary 
approach using Starburst99 models, in which the O star content is inferred from the 
stellar continuum, and the WR populati on is obtaine d from spectral synthesis of op- 
tical WR features using the grids from ISmith et all We find reasonable agreement 
between the two methods for the O star content and the N(WR)/N(0) ratio but find 
that the WR subtype distribution is in error in the Starburst99 models, with far too 
few WN stars being predicted. We attribute this failure to the neglect of rotational 
mixing in evolutionary models. Our various modelling approaches allow us to measure 
the cluster masses. We identify Al as a super star cluster (SSC) candidate with a mass 
of ~ 1-2 x 10 5 M . A total mass of ~ 6 x 10 5 M Q is inferred for the ionizing sources 
within Tol89-B. 

Key words: stars - Wolf-Rayet: stars - O stars: galaxies - starbursts: galaxies - 
massive star population: giant Hll regions - Tol 89. 
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1 INTRODUCTION 

Giant Hll regions (GHRs) are characterised by their 
large sizes (up to » 1 kpc), superso nic gas motions 
iMelnick. Tenorio-Tagle fc TerlevichlllflOflh and high Ha lu- 
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minosities (10 38 -10 41 ergs s ~ 1 : iKennicuttll 984^1 . A review of 
their properties is given riv lShieldTi^9(lfi . The presence of 
GHRs in galaxies denotes sites of recent, intense episodes 
of massive star formation. The nearest cxtragalactic GHR 
is 30 Doradus in the Large Magellanic Cloud (LMC) which 
hosts a compact sta r cluster (R136) of mass ~ 2-6 x 10 4 Mq 
llHunter et al.ll99fih as the main ionizing source. Conversely, 
another nearby GHR, NGC 604 in M33, is of similar size 
but contains multiple OB associations rather than a central 
massive cluster. The mechanism(s) which ultimately deter- 
mines whether a cluster or a complex of OB associations 



2 F. Sidoli et al. 



is for med when an intense star-for ming event occurs is un- 
clear. |Mme^ee^^^Efremo3 ((^23) suggest that clusters are 
formed in preference to loose associations in high pressure in- 
terstellar envir onments. Clusters wit h masses as high as 10 5 
to ~ 10 8 M l|Maraston et al.ll2004T) - often termed young 
massive clusters (YMCs), or super star clusters (SSCs) - are 
observed in the extreme en vironments of st arburst galax- 
ies and galaxy mergers (e.g. lWhitmorell2003h : the most lu- 
minous and compact of which show similar properties to 
older globular clust ers (GCs) such as those observed in 
the Milky Way (e.g. [ Holt zman et al lllQ92t IWhitmore et alJ 
1997USchweizer fc Seitzerll998l : lBastia,n et alJl2006l . see also 
Larsenl 2004 for a review). This has led to the suggestion 
that YM Cs may represent the young counterparts of old 
GCs (e.g. I Ashman fc Zer>dll992T) and may offer insight into 
the formation and evolution of GCs in the lo cal universe. 

Recently, IChen. Chu fc Johnson! ^OOfifl have studied 
the cluster content of three GHRs (NGC5461, NGC 5462 
and NGC 5471) in M 101. They find that they contain clus- 
ters similar to R136 in mass rather than the more massive 
SSCs, although they note that the three R136-like clusters 
in NGC 5461 may merge to form an SSC. Overall, they find 
evidence for a link between the molecular cloud distribu- 
tion and the cluster luminosity function in the sense that 
a diffuse distribution produces more lower mass clusters, in 
support of the hypothesis that massive clusters are formed 
in high pressure environments. 

The giant Hn region Tol89 l|Smith et al.lll976T ) is lo- 
cated at the south-western end of the bar in the late-typ e 
barred spiral (Sdm) galaxy NGC 5398 jPurret et alj fl9851. 
A Digital Sky Survey image of this galaxy is shown in Fig.^ 
Tol 89 is conspicuous in being the only large massive star- 
forming complex present in the entire galaxy. Tol 89 has an 
extent of « 24" x 18" or 1.7 x 1.2 kpc and an absolute blue 
magnitude of —14.8, (assuming a distance of 14.7 Mp c based 
on Ho = 75 km s^Mpc -1 ; ISchaerer et all 119991) which 
makes it one of the most impressive GHRs known. 

The presenc e of Wolf-Rayet (W R) stars in Tol 89 was 
first reported hv lDurret et al.Nl985ll who detected a broad 
emission bump at A4650. ISchaerer. Contini fc Kunthl (1999) 
identified broad features of N ill 4640, Hen 4686, and very 
strong C iv 5808, suggesting the presence of both late-type 
WN (WNL) and early-type WC (WCE) stars in Tol 89. From 
a detailed spatial analysis, they find an offset between the 
continuum and nebular lines, suggesting a complex star- 
forming region. I.Tohnson. Indebetouw fc Pisanol l|2003l) ob- 
tained radio observations of Tol 89 and discovered an unre- 
solved thermal radio source with a Lyman continuum flux 
of ~ 4500 x 10 49 s _1 , equivalent to an SSC with a mass of 
10 6 Mq if the source is a single cluster. In terms of radio lu- 
minosity, they find that Tol 89 is among the most luminous 
radio Hn regions so far observed and is comparable to the 
GHR NGC 5471 in M 101. Finally, Tol 89 w as one of the ob- 
jects in the survey of lChandar et al J l)2004h to measure the 
WR content of actively star-forming regions through ultravi- 
olet (UV) spectroscopy obtained with the Hubble Space Tele- 
scope (HST) Space Telescope Imaging Spectrograph (STIS). 

In this paper, we examine the massive stellar content 
of Tol 89 through an analysis of archival HST images and 
UV spectroscopy, and high resolution optical spectroscopy 
obtained with the Very Large Telescope (VLT). We show 
that Tol 89 is a young, very massive star-forming complex 




Figure 1. Digital Sky Survey R-band image of NGC 5398 show- 
ing the location of the GHR Tol 89. At the adopted distance of 
14.7 Mpc, 1' corresponds to ~ 4 kpc. The field of view is 5 x 5 
arc minutes. 



with at least seven young compact massive clusters. The 
fact that Tol 89 is located at the end of the bar in NGC 5398 
indicates it may have been formed t hrough gas inflow in a 
high pressure environment, although I.Tohnson et alJ l|2003tl 
suggest that the weak bars found in late-type galaxies are 
not strong enough to generate the required gas inflow. 

The paper is structured as follows. The observations 
and reductions are presented in Section [2] In Section |^1 we 
describe the spectra and in Sections 0] and we derive the 
properties of the knots and their ionizing clusters. In Sec- 
tion HJ we estimate the massive star content using empirical 
and synthesis techniques. Finally, in Sections and El we 
discuss and summarise our results. 



2 OBSERVATIONS 

We have obtained high spectral resolution UV- Visual 
Echelle Spectrograph (UVES)+VLT echelle spectroscopy of 
the brightest optical knots in the giant Hn region Tol 89. 
We supplement this dataset with archival HST imaging 
and spectroscopy obtained using the STIS CCD and FUV- 
MAMA detectors in the optical and UV respectively (Pro- 
posal ID 7513; C. Leither er, P.I.). The UV spe ctroscopic 
data are also presented in IChandar" et al.l l)2004l) . A sum- 
mary of the observations is given in Table U All the HST 
archive data were processed using the standard CALSTis data 
reduction pipeline. 

2.1 Imaging 

2.1.1 Reduction 

Optical and UV images were obtained on 1999 February 19 
using the MIRVIS/STIS CCD and MIRFUV/FUV-MAMA 
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Figure 2. STIS MIRVIS (A con =7230 A) image showing the po- 
sition and orientation (PA = -90°) of the optical STIS slit, of size 
52" x O'.'l arc seconds. Clusters Al-3, Bl, CI and C2 are indicated 
(see Section \2. 1.21 . At the adopted distance of 14.7 Mpc, 1" cor- 
responds to ~ 70 pc. The stretch of the image has been set so as 
to highlight the position of cluster A3. The insets show the regions 
around Al and A2 (top right) and Bl (centre left) at different 
stretches to highlight better these clusters. 




Figure 4. STIS MIRVIS (A con =7230 A) image showing the posi- 
tion and orientation (PA = 90°) of the UVES slit, of size l'.'4x 10". 
Knots A, B and C are indicated (see Section l2.1.2l . The orientation 
and scale are as in Figs. l2l and l3l 



longpass filters/detectors. The optical dataset consists of 
four image sets. Before combining, we first used intensity 
histograms of the images to determine a threshold value for 
the cold pixels, below which all pixels were flipped to high 
values. These were subsequently removed, along with cos- 
mic rays and hot pixels, using the iraf/stsdas cosmic ray 
rejection (occrej) task upon combining the images. The re- 
sulting images were then geometrically corrected using the 
x2d task. Remaining hot pixels were removed using the 
Laplacian cosmic r ay identification algorithm, L. A. Cosmic 
l|van Dokkumll200lf) . In Figs. |2| and El we show the optical 
and UV STIS images of Tol 89, with the slit positions of the 
STIS spectroscopy superimposed (see Section Pi. 2. 111 . Simi- 
larly in Fig. 131 we show the optical STIS image with the slit 
position of the UVES spectroscopy overlaid (Section l2.2.2ll . 

2.1.2 Optical and UV morphology 

Tol 89 is a young massive star forming complex as shown 
by the optical and UV STIS images in Figs. PJ and PJ It 




Figure 3. STIS MIRFUV (A ccn =1480 A) image showing the po- 
sition and orientation (PA = -90°) of the UV STIS slit, of size 
52" x Clusters Al-4, Bl, CI and C2 are indicated (see Sec- 
tion l2.1.2l . At the adopted distance of 14.7 Mpc, 1" corresponds to 
~ 70 pc. The stretch of the image has been set so as to highlight 
the position of clusters A3 and A4. The inset shows the region 
around clusters Al and A2. Note, cluster Bl is almost completely 
obscured in the UV. 

comprises three distinct knots of star formation that we label 
A, B and C in order of decreasing (optical) brightness. With 
the high spatial resolution of HST, we are able to resolve 
each individual knot into multiple cluster components. 

Along the optical STIS slit, the boundary between A 
and B is defined by the mid-point between the brightest op- 
tical component in each knot. We then identify the clusters 
in order of decreasing brightness along the slit, with 1 being 
the brightest (see Figs. Pil and PJJl. The UV counterparts to 
each cluster were then identified in the FUV-MAMA image. 
The separation between the clusters, as measured from the 
optical STIS image, are as follows: A4-*A1: l'.'53 (107 pc); 
A1-»A2: 0'.'32 (22 pc); A2^A3: 0'.'69 (49 pc); and A3^B1: 
l'/ 5 1 (106 pc). The projected separations in parsecs are given 
in parentheses, where 1" ~ 70 pc (assumin g a distance of 
14.7 M pc based on H = 75 km s~ 1 Mdc~ 1 : ISchaerer et alJ 
Il999h . 

Knot A is made up of multiple clusters that are bright 
in both the optical and the UV, and dominated by a young 
compact cluster which we denote as Al. Whilst knot B is 
also optically bright and complex it is much fainter than A 
in the UV. It is also the lo cation of the impress ive thermal 
radio source discovered bv l.Tohnson et ail l)2003h . From the 
STIS UV image (Fig. PJJ we can see that knot C contains 
two young compact massive clusters which we denote as CI 
and C2; where CI is the location of th e north-eastern spur 
of thermal radio emission identified bv l.Tohnson et alJl2003l 
(see their fig. 4). 

iDurret et all identify two regions in Tol 89 that 

they denote as X and C, where C is their centre of Tol 89. We 
identify our knot C with their region X and their centre of 
Tol 89 with our knots A and B. Within the centre of Tol 89, 
ISchaerer et alJ l|l999f) identify an area with maximum neb- 
ular intensity ( = our B), offset by ~ 2 arc seconds to the 
east of the region with maximum continuum intensity and 
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Table 1. Summary of the observations. The HST+STIS imaging and spectroscopy (Proposal ID 7513, C. Leitherer, P.I.) have been 
obtained from the HST archive. 
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HST+STIS Imaging, 1999 February 19 
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a Resolutions and dispersions are for CCD1/CCD2 



the location of the WR emission ( = our A; see also their 
fig. 5 1 ). 



2.2 Spectroscopy 

2.2.1 STIS spectroscopy 

The two-dimensional STIS spectroscopic data were obtained 
over two visits on 2000 March 20-21 using the G430L and 
G750M gratings and the STIS CCD detector for the op- 
tical regime, and the G140L grating and the STIS FUV- 
MAMA detector for the UV. The slits were centred on the 
brightest cluster (=A1) at co-ordinates of a = 14 ,1 01 m 19 f :92; 
S = -33 h 04 m 10 , :7 (J2000). The respective slit sizes are 
52" x O'.'l and 52" x 075 and are shown in Figs. El and 
Further details are given in Table U The image sets (two 
for each optical grating and four for the UV) were com- 
bined using the method described in Section 12.1.1 1 and rec- 
tified, wavelength and absolutely flux calibrated using the 
x2d task. From measurements of the FWHM of Gaussian 
fits to the unresolved H/3 and C n 1335 lines in the G430L 
and G140L spectra of Bl, we obtain spectral resolutions of 
4.9A and 3.lA respectively, or 1.8 pixels, giving a resolution 
of I.OA for the G750M grating. 



1 Note that the orientation in fig. 5 of lSchaerer et al] 1999 is in- 
correct. North points towards the top right corner of their middle 
panel rather than the top left. Thus, the peak in nebular inten- 
sity (our knot B) is offset by ~ 2 arc seconds to the east of the 
continuum intensity peak (our knot A) rather than south, as is 
suggested in their figure. 



From the UV spectral image, we identify five regions 
for extraction which we denote as Al, A2, A3, A4 and Bl in 
Fig- El (see Section 12.1.2(1 . The extraction widths correspond 
to 0.23, 0.12, 0.19, 0.19 and 0.30 arc seconds, respectively. 
Note that in Fig.EJthe bright UV source in the vicinity of Bl 
is actually the UV counterpart to the optical point source 
that lies below, and just outside of the narrower optical slit 
(see leftmost inset of Fig. EJ- For this reason we make care- 
ful considerations regarding the extraction width of Bl to 
ensure that there is as little contamination from this source 
as possible. 

For the optical we extract three regions Al+2, A3 and 
Bl, with the extraction widths being defined by the full 
width at ~10 per cent of the peak of the intensity profile 
at Ha. These extraction widths, corresponding to 0.55, 0.40 
and 1.13 arc seconds respectively, were then applied to the 
G430L grating spectra. Note, regions Al and A2 are barely 
resolved in the optical STIS 2D spectral image and hence are 
extracted as a single source object while cluster A4 is not 
detected above the background noise. Background subtrac- 
tion was performed by selecting regions of sky free of nebular 
emission. No correction for slit losses have been made. 



2.2.2 UVES spectroscopy 

Echelle spectra of Tol 89 were obtained in serv ice mode on 
2004 May 6 with UVES JD'Odorico et alJl2000h at the VLT 
Kueyen Telescope (UT2) in Chile (Proposal ID 73.B-0238A, 
L. J. Smith, P.I.). UVES is a two arm cross-dispersed echelle 
spectrograph with the red arm containing a mosaic of an 
EEV and a MIT-LL CCD. The blue arm has a single EEV 
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CCD; all three CCDs have a pixel size of 15/xm. Simultane- 
ous observations in the blue and the red were made using 
the standard setups with dichroic #1 (346+564 nm) and 
dichroic #2 (437+860 nm), covering an almost continuous 
wavelength region from 3100-10 360A; the regions between 
5610-5670A and 8540-8650A were not observed as a result 
of the gap between the two CCDs in the red arm. The slit 
was oriented at a position angle of 90 deg to pass through 
the two brightest knots in Tol 89, hereafter A and B when 
discussing the UVES data (see Fig. EJ. The slit width was 
set at 1.4, giving a resolving power of ~ 30 000 in the blue 
and ~ 28 000 in the red. The pixel scales for lxl binning 
are shown in Tabled 

At the time of the observations the conditions were clear 
and the seeing was typically better than l'/2, A total inte- 
gration time of 47.5 min. was divided into two equal length 
exposures to prevent saturation of the brightest nebular lines 
(Ha, H/3 and the [O ill] AA4959, 5007). Observations of the 
spectrophotometric standard Feige 67 were also made using 
the same set up for the purpose of performing flux calibra- 
tions. Spectra of Th-Ar comparison arcs were obtained to 
perform wavelength calibrations. 

The data were reduced according to the standard CCD 
and echelle reduction procedures within iraf using the cc- 
dred and echelle packages. This included the subtrac- 
tion of a bias level determined from the overscan region, 
bias subtraction, division by a normalised flat field and bad 
pixel correcti on. Images were cle aned of cosmic rays using 
L. A. Cosmic l|va,n Dokkumll200lh . The echelle orders were 
extracted using widths set to 2.4 and 3.6 arc seconds for 
knots A and B respectively and wavelength calibrated, at- 
mospheric extinction corrected and flux calibrated. Resolu- 
tions - as measured from the FWHM of Gaussian profile fits 
to the Th-Ar arc lines - and dispersions are given in Tabled 
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Figure 5. Velocity corrected STIS FUV-MAMA G140L spec- 
tra showing the extractions of cluster Al, A2, A3, A4 and Bl 
(top, 2nd, 3rd, 4th and bottom panels respectively). P Cygni pro- 
files due to the winds of massive O and WN stars are identified: 
Si iv 1400 (O supergiants), Nvl240 and Civ 1550 (O stars) and 
He n 1640 (WN). In the case of A2, the strength of the C iv 1550 
feature, which exceeds that of He n 1640, suggests the presence of 
WC stars. 



3 DESCRIPTION OF THE SPECTRA 
3.1 STIS UV 

The velocity corrected (V hc i w 1230-1240 km s" 1 ; see Sec- 
tion UV spectra of clusters Al-4 and Bl are shown in 
Fig. the characteristics of stars with strong stellar winds 
due to the presence of P Cygni profiles of Si iv 1400, N v 1240 
and C iv 1550 can be seen. In the case of A3, the narrow and 
resolved Si iv 1400 profile is suggestive of late-type O super- 
giants. 

Since O stars do not show strong He n 1640 emission, 
its presence in Al, A2 and A4 can be attributed to the winds 
from WN stars. In A2 the emission strength of the C iv 1550 
feature, which is greater than He n 1640, suggests that WC 
stars are also present, while the non detection of Hen 1640 
emission in A3 and Bl would suggest an absence of WRs. 
In cluster A4 we appear to detect N iv] 1486 emission from 
mid-WN stars and note that its UV spectral appearance 
bares a close resemblance to that of a mid-WN star (see 
fig. 2 oflCrowther fc DessartllT99l) . 

The UV S TIS data are also presented in fig. 1 of 
IChandar et alJ l)2004l) . who used a larger aperture width of 
81 pixels to obtain their ID spectrum of the region they de- 
note as Tol 89-1. Their extraction comprises our regions Al 
and A2, plus diffuse inter-cluster emission and also shows 



the P Cygni profiles of Nvl240, Si iv 1400, and Civ 1550, 
as well as He II 1640 emission from WN stars. Our smaller, 
individual extractions of Al and A2 (~ 9 and 5 pixels re- 
spectively) have allowed for different WR populations (WN 
and WC) to be identified in the individual clusters. From 
the spectrum of Tol 89-1 shown in fig. 1 of IChandar et alJ 
the presence of WC stars is ambiguous and thus only a WN 
population can be inferred from their extraction. 



3.2 STIS optical 

The optical spectrum of Al+2 shows a very blue continuum 
with very weak emission lines of Ha, H/9, [O m] 4959, 5007 
and [N n] 6548/84. The portion of the G430L spectrum in the 
region of He II 4686 is shown in Fig.UJi; we identify a possible 
nebular Hen 4686 component. The FWHM is comparable to 
that of the unresolved [O m] 5007 emission line (cf 7.0+1.8 
and 6.7+1.0 A respectively). We briefly discuss the origin of 
this emission in Section 

Clusters A3 and Bl also show emission lines of Ha, 
H/3 and [Om]4959, 5007; although they are weak in A3. 
In addition to these lines, Bl shows weak emission lines of 
[O n] 3727, Hy and H<5. Both A3 and Bl are significantly red- 
dened compared to Al+2; no WR emission features are de- 
tected, in agreement with the UV STIS spectra (see Fig. . 
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Figure 6. Velocity corrected STIS G430L spectrum of Al+2 
(top) and UVES+VLT spectrum of B (bottom) showing possible 
nebular He n 4686 emission. Nebular lines of H/3 and [Oin] in 
A1+A2, and He i , [Fern] and [Anv] in region B are also identi- 
fied. 



3.3 UVES 

The high S/N, high spectral resolution UVES data for re- 
gions A and B are shown in Figs. and Irrespectively over 
the wavelength range of 3600-7400 A. The spectra are rich 
in emission line features typical of Hn regions. The strongest 
emission lines, those of Ha, H/3, [O n] 3726/29, 7319/30, 
[O m] 4959, 5007, [N ii] 6548/84, and [S n] 6716/31, exhibit 
a broad component resulting in extended wings. In the case 
of region A, the narrow components of these emission lines 
are split into two velocity components (denoted as VI and 
V2). Both spectra show little or no stellar absorption in the 
Balmer lines. Metal lines due to an older population of stars 
are also absent and imply that regions A and B consist of a 
young population of stars with very similar ages. An analysis 
of the nebular lines is presented in Section 

Reg arding stellar fe atures, we confirm the observations 
made by ISchaerer et all detecting broad WR emissions in 
knot A at A4640 (the blue bump) and A5808 (the yellow 
bump), see Figs. and |2| The blue bump is a blend of neb- 
ular [Fe m] 4658 and stellar Nv4620, Nm4640, Cm/Civ 
A4650/58 and Hen 4686, suggesting the presence of both 
early (WNE) and late-type (WNL) nitrogen-rich WR stars. 
The equivalent width of the blue bump is ft! 11 A, with A4686 
contributing « 3A. We also detect broad (FWHM ~ 80 A) 
C iv 5808 emission from WC stars with an equivalent width 
of ~ 12 A. 

In region B only nebular lines of [Fem]4658, 4701, 
[Ar iv] 4711/40 and He 1 4713 are present in the region of 
the blue bump, see Fig.OJj. No broad WR emission is seen. 



Again, we detect the presence of nebular Hen 4686. The 
width of He II 4686 emission is comparable to that of [Fein] 
and [Anv] , cf 0.6 ± 0.2 and 0.7 ± O.lA, respectively. The 
origin of this nebular He n 4686 emission is discussed in 
Section 



4 PROPERTIES OF THE KNOTS AND THEIR 
IONIZING CLUSTERS 

In this section we derive the ages and extinction towards 
knots A and B and their ionizing clusters Al-4 and Bl using 
stellar (Section 14.111 and nebular (Section 14.211 diagnostics 
in the UV and optical respectively. Cluster sizes for Al-4, 
Bl, CI and C2 are also derived from the UV STIS image 
in Section 14.1.11 Note that no age or extinction estimates 
can be made for CI and C2 since no spectroscopic data is 
available for these clusters. 



4.1 UV derived cluster properties 

4-1.1 Sizes 

Cluster sizes were det ermined from t he STIS UV image using 
the ishape routine of lLarsenl l|l999f) . The routine compares 
the observed cluster profile with PSF-convolved analytical 
functions that model the surface brightness of the clusters. 
The best-fitting function, determined by the % 2 results, then 
gives the effective radii R a s, or half light radius for each 
cluster. A correction for the true effective radii for elliptical 
profiles -Rcff.cii, as given by equation 11 of the ishape user's 
guide, is made. 

We choose to model our clusters using lKingl jl962T ) pro- 
files for different concentration parameters c - where c is 
defined as the ratio of the tidal radius r t to the core radius 
r c and takes values of 15, 30 and 100 - and moffat pro- 
files for different power indices a, where a equals 1.5 or 2.5 
(moffat15 or moffat25). Various fitting radii were used, 
ranging from a minimum radius that is set equal to the ap- 
proximate size of the cluster to a maximum radius that is 
set to be the furthest one can move from the cluster with- 
out background contamination from a nearby source. The 
results of the models giving the most internally consistent 
FWHM and best \ 2 statistics are shown in Table El All the 
clusters appear to be very co mpact, with h alf light radii < 3 
pc and sizes typical of SSCs l)Larsenll200l . 



4.1.2 



We a dopt the method of IChandar. Leitherer fc Tremontil 
i2004l) for determining the ages of our clusters by first 
correcting the STIS UV spectra for G alactic foreground 
extinction. We use the ISeatonl lll97<t) Galactic law for 
an_E(B-F) ga i = 0.066 nw JSchWel. Finkbeiner fc David 

1998) . We then normalise the spectra and c ompare them to 
the m odel data output of Starburst99 v5.0 l|Leitherer et alJ 

1999) sc aled from a 10 6 M Q instantaneous burst with a 0.1- 
100 Mq iKroupal l)200l|) IMF. 

Based on our abundance determinations made in Sec- 
tion 15.21 an LMC metallicity was assumed. UV O star 
spectra shortward of A < 1600 A are taken from tem- 
plate stars in the Large and Small Magellanic Clouds 
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Figure 7. Velocity corrected UVES+VLT spectra of Tol 89-A covering the spectral range of 3600-7400 A. In each panel the spectra are 
scaled to arbitrary flux values to show the detail in the spectra. The unsealed spectra are plotted to show the relative line intensities of 
the strongest emission lines. The velocity splitting in the strongest emission lines can be seen. The blue and yellow WR bumps are also 
indicated. 
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Figure 8. Velocity corrected UVES+VLT spectra of Tol89-B covering the spectral range of 3600-7400 A. In each panel the spectra are 
scaled to arbitrary flux values to show the detail in the spectra. The unsealed spectra are plotted to show the relative line intensities of 
the strongest emission lines. 
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Table 2. Results of the best-fitting analytical profiles giving the most internally consistent x 2 statistics and FWHM for clusters Al-4, 
CI and C2 seen in the UV STIS image. The FWHM and minor-to-major axis ratios are averages of between 4-6 measurements over 
different fitting radii. 



Cluster 


Model 


FWHM 


minor 


^eff 








(pix) 


major 


(pix) 


(pix) 


(") 


(PC) 


Al 


MOFFAT25 


2.4 


0.79 


1.63 


1.46 


0.036 


2.6 


A2 


Kingl5 


2.4 


0.73 


1.69 


1.46 


0.036 


2.6 


A3 


MOFFAT25 


2.8 


0.75 


1.92 


1.68 


0.041 


3.0 


A4 


Kingl5 


1.0 


0.76 


1.05 


0.92 


0.023 


1.6 


CI 


MOFFAT25 


2.5 


0.51 


1.67 


1.26 


0.031 


2.2 


C2 


MOFFAT25 


2.0 


0.71 


1.37 


1.17 


0.029 


2.1 





u 

CIV 




UVES 

Tol89-A 'yellow 1 WR bump 



Wavelength (A) 

Figure 9. Velocity corrected UVES+VLT spectra showing the 
blue (top) and yellow (bottom) WR bumps in Tol89-A. The 
blue bump is a blend of nebular [Fem]4658 and stellar Nv4620 
Nm4640, Cm/Civ A4650/58 and Hen4686 - indicating the 
presence of both early (WNE) and late-type (WNL) nitrogen rich 
Wolf-Rayet stars. The yellow bump is due to C iv emission from 
the strong stellar winds of WC4-5 stars. 



(Le itherer et "ai"l l200ll) . All other empirical data are from 
spectral types at solar met a llicity , i.e. the atlases of 
Robert. Leitherer fc Heckmanl lll993f) for WR stars and 
de Mello. Leitherer fc Heckmar]T20001) for B stars. 

In Fig. ED we show the best-fitting models for clusters 
Al and A2. We can see in the case of Al that we are able to 
reproduce well the strengths of the Si iv 1400 and C iv 1550 
O-star P-Cygni profiles, enabling us to constrain the age to 
4.5±0.5 Myr. For A2 we determine a similar age of 5±1 Myr. 
The uncertainty is larger since we rely solely on the fit to 
Si iv 1400. This is due to the fact that the Starburst99 mod- 
els fail to reproduce the C iv 1550 feature, whose strength 
implies a significant contribution from WC stars in A2; un- 
like in Al where the fit to Civ 1550 is excellent. The WR 
spectral synthesis in the UV is rudimentary because it relies 



on solar metallicity spectra of a few WR stars. Until detailed 
WR spectral synthesis is available for the UV wavelength 
range, it is difficult to draw any quantitative conclusions 
on the WR content of regions Al and A2 from Starburst99 
modelling alone. Our derived ages for Al and A2 are consis - 
tent with the 4±1 Myr determined bv lChandar et al.1 <2004h 
for their extraction of Tol89-l (= our Al + A2). 

Our age estimates for A3, A4 and Bl are somewhat 
more uncertain given the poorer S /N of the spectra; we ob- 
tain approximate ages of 3.0-5.5 Myr for A3 and A4 and 
< 3 Myr for Bl from the model fits. The age estimate for 
A4 is especially uncertain, due to the weakness of Si iv , with 
mid-WN stars contributing to the C iv profile, such that 
Starburst99 modelling is inadequate. We suspect that A4 
has an age of ~3 Myr, given its similarity to mid-WN stars 
observed in young massive clusters f Section Ifi. 1.2(1 . 

Given the age of cluster A3 one would also expect WR 
stars to be present. Whilst WC stars do not appear to be 
present from the observed CIV 1550 profile, weak Hell 1640 
from WN stars is not excluded given the observed low S/N 
for this region. The presence of WN stars in A4 is discussed 
in Section 



4-1.3 Extinction 

From the best-fitting age models we determine an estimate 
for the internal reddening E(B-V)i n t towards each cluster 
by dereddening the Galactic foreground extinction-corrected 
cluster spectra to match the slope of the best fitting Star- 
burst99 model. The fit is performed over the wavelength 
range of 1240-1600 A and the models normalised to match 
the continuum either side of C iv 1550. We adopted the stel 
lar LMC extinction law of[ 



owa: 



to der edden our 

UV spectra in favour of the C,alzett^t ajj l|20 0(i[) starburst 
obscuration law used bv lChandar^tHjil l)2004|) . The latter is 
more appropriate for unresolved star forming galaxies than 
for effective point sources, such as is the case for Tol 89-A1 
an d A2. Similar conclusi o ns we re reached for NGC3125-A1 
bvlHadfield fc Crowtherl EoM). 

The following internal reddening values are obtained 
for Al and A2 respectively: 0.09 ± 0.02 and 0.08 ± 0.02 
mag. This is equivalent to ~ 1.0 ± 0.2 and 0.9 ± 0.2 mags 
of int ernal extinction at A1500 (^4isoo) for the Irlowarthl 
lll983h LMC reddening law. Although our reddening esti- 
m ates are in agree ment with the value of 0.08 mag obtained 
bv lChandar et all the extinction at A1550 o btained here is 
approximately twice that obtained using the lCalzetti et all 
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Figure 10. The rectified FUV spectra (solid black line) of clus- 
ters Al (top) and A2 (bottom) shown with the best-fitting (age) 
Starburst99, LMC models (dotted line). The STIS spectra have 
been smoothed for clarity and velocity corrected. 



starburst obscuration law; for which A1500 = 0.4 mag for 
E(B-V)int— 0.08 mag. For regions A3, A4 and Bl we obtain 
E(B- V)int of 0.10 ± 0.03, 0.08 ± 0.03 and 0.17 ± 0.03 mag, 
or A1500 ~ 1-1 ± 0.3, 0.9 ± 0.3 and 1 .8 ± 0.3 mag respec- 
tively. Again, using the lCalzetti et all law the extinction at 
A150 is approxima tely half that obtained using the LMC 
law jHowartblll983h . 



4.2 Optically derived knot properties 

4.2.1 Ages 

Estimates for the ages of knots A and B were determined 
from the UVES spectra by comparing the equivalent widths 
of nebular H/3 emission to Starburst99 predictions. We 
obtain equivalent widths of 55 and 295 A for A and B 
respectively. For an instantaneous burst of LMC metal- 
licity and a 0.1-100 Mq Kroupa IMF, this corresponds 
to ages of ~4.5 and 2.5 Myr; assuming that all hydro- 
gen ionising photons are absorbe d within the Hn region. 
ISchaerer. Contini fc Kunthl l|l999fl determine an age of 4.5- 
5.0 Myr for Tol 89 by comparing their measured H/3 equiva- 
lent widths with standard, SMC population synthesis mod- 
els. 

The ages of the knots derived from the EW(H8) are in 
excellent agreement with the cluster ages derived in Sec- 
tion 14.1.21 The average cluster age in knot A (Al-4) is 
~ 4.5 ± 1.0 Myr, which is in excellent agreement with the 
4.5 Myr determined for Tol 89-A. An age < 3 Myr is derived 
for Bl, which is also in good agreement with the < 2.5 Myr 
derived for Tol89-B. 



4-2.2 Extinction 

In the optical, the total extinction E(5-U)tot 2 is determined 
from the Balmer line decrement. For both the STIS and 
UVES data we use only the Ha and H/3 line fluxes to deter- 
mine the extinction for the following reasons. In the STIS 
data, the higher order Balmer lines of H7 and H<5 are very 
weak and thus are not reliable indicators of E(B-V) to t. In 
the UVES data we observe little to no underlying stellar ab- 
sorption. After correction for galactic foreground extinction 
we determine an internal reddening value for knot A which is 
consistent with there being zero reddening. The low extinc- 
tion towards knot A is surprising given it's young age (~ 4.5 
Myr) and suggests that clusters very quickly disperse their 
natal clouds within a few Myr. For knot B, we obtain an 
internal reddening value E(B- F)int of ~ 0.24 from both the 
STIS and the VLT spectra. We therefore adopt the follow- 
ing (total) reddening values for knots A and B of 0.07±0.01 
and 0.29 ± 0.03 mag respectively. Previous determinations 
of the average reddening over the Tol 89 complex lie within 
the bounds of our estima tes, cf 0.12 jTerleyich et alJll99ll : 
C{HB) = 0.18) and 0.20 jPurret et al.lll985h . 



5 NEBULAR PROPERTIES OF THE KNOTS 

Line fluxes and equivalent widths for nebular lines were de- 
termined using the elf (emission line fitting) and ew (equiv- 
alent width) routines within the starlink package dipso. 
Fluxes were measured by Gaussian fitting, allowing for line 
centres and widths to vary freely - except when fitting dou- 
blets, in which case widths and relative line centres were con- 
strained. The observed and intrinsic line fluxes normalised 
to H/3 = 100 are listed in Table El 

We find that the strongest nebular emission lines are 
composed of a narrow and a broad component, both centred 
at similar velocities. In addition for knot A, we detect line 
splitting in the narrow component with a separation between 
the blue (VI) and red (V2) components of 48 ± 1 km s _1 . 
The FWHMs of the broad and narrow components are for 
knot A: 111 ±2, 32 ± 2 (blue component), 41 ± 2 km s" 1 
(red component), and 71 ±2, 27 ±2 km s _1 for knot B. The 
multiple components are seen in the following emission lines: 
Ha, H/3, [On] AA 3726/29, AA 7319/30, [O m] AA 4959, 5007, 
[No] AA 6548/84, and [S 11] AA 6717/31. In lines with poorer 
S/N no broad or velocity-split components are detected. 

In Fig. llll we show, as an example, the [S 11] AA 6716/31 
emission seen in A (top) and B (bottom) with the best-fitting 
Gaussians superimposed. The broad components account for 
a significant fraction of the total line flux (31 and 50 per cent 
for A and B respectively). The broad components are sym- 
metrical with respect to the mean velocities of the narrow 
components and have the same central velocities for both 
A and B, suggesting a common origin. The only apparent 
difference between the broad components in the A and B 
spectra is in the line width (111 vs. 71 km s _1 ). 

Mean heliocentric radial velocities for A and B were 
determined by averaging the velocities of all the Gaussian 

2 Here we define E(B-V)tot to be the Galactic foreground con- 
tribution E(B-V) gs .i (0.066 mag; Schlegel et al. 1998) plus the 
internal contribution, ~E(B-V)i nt . 
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Table 3. Observed (Fx) and intrinsic (Ia) nebular line fluxes 
for knots A and B relative to H/3 = 100, and equivalent widths, 
Wj, of WR emission lines seen in knot A. No WR emissions are 
detected in knot B. 



Ao 


Ion 


Fa 


Ia 


Fa 


Ia 


Nebular Emission Lines 


3726 


[On] 


73.92 


78.23 


39.34 


50.46 


3729 


[On] 


101.20 


107.10 


51.43 


65.93 


3869 


[Nem] 


26.33 


27.70 


27.95 


34.95 


3967 


[Nem] 


6.63 


6.94 


7.25 


8.90 


4026 


He i 


1.81 


1.90 


1.50 


1.81 


4069 


[Sn] 


0.74 


0.77 


0.71 


0.85 


4076 


[Sn] 






0.18 


0.21 


4102 


H<5 


22.49 


23.41 


22.08 


26.34 


4121 


He i 






0.18 


0.22 


4144 


He i 






0.18 


0.21 


4341 


H7 


48.87 


50.25 


41.84 


47.31 


4363 


[Om] 


2.28 


2.34 


2.60 


2.92 


4471 


He i 


3.22 


3.29 


3.40 


3.73 


4658 


[Fein] 


0.68 


0.68 


0.57 


0.60 


4686 


He ii 






0.12 


0.12 


4702 


[Fein] 






0.10 


0.10 


4711 


[Ar iv] 


0.50 


0.51 


0.22 


0.23 


4713 


He i 






0.45 


0.46 


4740 


[Ar iv] 


0.32 


0.33 


0.24 


0.25 


4959 


[Om] 


121.40 


120.80 


166.40 


162.50 


5007 


[Om] 


363.20 
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498.50 
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5016 


He i 


2.18 


2.16 


2.06 


1.99 
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He i 






0.11 


0.10 


5518 


[Clin] 


0.52 


0.50 


0.49 


0.42 


5538 


[CI m] 


0.41 


0.39 


0.31 


0.26 


5755 


[Nil] 






0.46 


0.38 


5876 


He i 


11.60 


11.07 


15.07 


12.24 


6300 


[Or] 


1.85 


1.74 


1.71 


1.30 


6312 


[Sm] 


1.42 


1.33 


2.02 


1.54 


6364 


[Oi] 


0.22 


0.21 


0.53 


0.40 


6548 


[Nii] 


4.53 


4.23 


5.42 


3.98 


6563 


Ha 


290.20 


270.40 


390.50 


286.30 


6583 


[Nil] 


14.91 


13.88 


17.26 


12.62 


6678 


He i 


3.09 


2.87 


4.80 


3.47 


6716 


[Sn] 


13.29 


12.33 


16.31 


11.71 


6731 


[Sn] 


9.97 


9.24 


12.80 


9.18 


7065 


He i 


2.44 


2.24 


4.29 


2.95 


7136 


[Arm] 


8.99 


8.24 


16.24 


11.06 


7281 


He i 


0.52 


0.47 


0.91 


0.61 


7319 


[On] 


1.56 


1.42 


0.79 


0.52 


7331 


[On] 


0.84 


0.76 


1.19 


0.79 


9069 


[Sm] 


27.70 


24.35 


56.23 


31.90 


9531 


[Sin] 


97.78 


85.32 


199.60 


109.70 


log H/3 




-13.32 


-13.23 


-12.86 


-12.43 


B(B-V)tot 




0.07 


0.29 
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Figure 11. UVES spectra of knots A (top) and B (bottom) 
in the spectral range of the [S n] 6716, 6731 doublet showing 
the multiple Gaussian component fits to the line profiles. The 
broad (long-dash), narrow (short-dash) and summed components 
(dotted line) are plotted. In region A, the two narrow velocity 
components, VI (short- dashed) and V2 (dash-dot-dot-dot line) 
can be seen. The central velocities of the Gaussian fits are given 
in each plot. 



components (broad and narrow) fit to the emission lines in 
each region and are found to be 1241±2 and 1234±2 km s _1 
respectively. These results are consi stent with the valu es re- 
ported in the liter ature: 1232 ± 51 fl Pu rret et alJll98fJl and 
1226 ± 11 km s" 1 ijSchaerer et alJll99<tl . 

Overall, we find that knots A and B have similar veloc- 
ities, differing by 7 km s _1 , and thus are part of the same 
star-forming complex. For A, we see velocity splitting, sug- 
gesting that the winds from the clusters have swept up the 
surrounding interstellar gas into a shell which is currently ex- 
panding at 24 km s _1 . For B, we see no line splitting which 
is in accord with its younger age. We also detect a broad 
component which appears to be common to both regions A 
and B. The maximum width of these features corresponding 
to the full width at zero intensity (FWZI) is 1140 and 940 
km s _1 for A and B; this suggests that the two knots contain 
gas with velocities up to ~ 450-600 km s _1 . We consider the 
origin of this high velocity component in Sect. 



5.1 Electron densities and temperatures 

Electron densities N c , and temperatures T e , were deter- 
mined using the diagnostic line ratios listed in Table 0] 
The individual broad and narrow velocity components were 
summed in our calculation of N e and T c due to the fact 
that no broad or velocity shifted components were detected 
for the weaker, poorer S/N, temperature-sensitive diagnostic 
lines (AA4069, 4076, 4363 and 5755). We note that the ratios 
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of the line fluxes for the broad and narrow components in 
[S n] are identical (see Fig. lilt suggesting that they have 
similar electron densities. 

Electron densities were calculated assuming a constant 
T c of 10 000 K. Since [O n] and [S n] give consistent N e (see 
TableEJ we take the average of the two to be representative 
of the densities in Tol 89-A and B. We find N c of « 90 ± 40 
and ~ 150 ± 40 cm -3 respectively. [Ar iv] and [CI in] probe 
denser regions of gas and thus give higher N e estimates. 
Refined estimates for the electron temperatures were then 
determined using the average N c determined from + and 
S + above. Using the most consistent T c between the different 
diagnostic line ratios, we derive the following 'average' T c of 
« 10 000 ± 300 K for A; neglecting T(S + ) due to its large 
uncertainty. 

For region B we derive a value of 9 800 ± 300 K. We 
exclude the large T c derived from N + and + because 
in moderate to high-excitation Hn regions the [N n] 5755 
and [O n] 7330 lines can be excited by the reco mbination of 
N 2+ and Q 2+ i n the higher-excitations zones 
iLiu et aliEoorl . thus leading to higher electron tempera- 
tures (see Table EJ. 



5.2 Abundances 

Ionic abundances were calculated using the average N e and 
T c determined in Section l5~Tl fsee Table HJ. To determine 
total abundances we adopt the following expressions given 
by Eqn.Q where the ionisation correction factors (ICF) for 
N, S and Ar are as defined in equations 15, 18 and 19, re- 
spectively, of llzotov. Thuan fc Lipovetskvl l|l994h . 



H 



o + + o 2 



- = ICF(S) 



S+ + s 2+ 

H+ 



N N + 



Ar 

— = ICF(Ar) 



Ar 2 + + Ar 3 + 
H+ 



(1) 



We derive 12+log(0/H) for knots A and B respectively 
of 8.27 and 8.38 w hich are in good agree ment with the value 
of^ derived bvlSchaerer etTbl B: whilelDurret et~aTl 
( 1985) derive a much lower oxygen content of 8.03. Oxy- 
gen and s ulphur abundances are very similar to those of 
the LMC ijR.ussell fc Donitalll99oT l. such that hereafter we 
assume an LMC metallicity for Tol 89, although nitrogen 
and argon abundances are somewhat lower. With respect 
to 30 Dor l|Peimbertll2003^ . knots A and B are moderately 
depleted by 0.25-0.5 dex. 



6 MASSIVE STAR POPULATION 

In the following section we examine the massive star popu- 
lation in Tol 89, applying both empirical and synthesis tech- 
niques in the optical and the UV. 



Table 4. Derived physical properties and abundances from the 
UVES spectra of Tol 89 for knots A and B. A c omparison is made 
with the abunda nces derived fo r the LMC iRussell Kr. Donital 
Il990l) and 30 Dor JPeimberfeOolh . 



Diagnostic 


Knot 




Line Ratio 


A B 


LMC° 30 Dor 6 



Density 

[On] A3726/A3729 
[S n] A6731/A6716 
[Ar iv] A4740/A4711 
[CI m] A5537/A5517 

Temperature 
[Nil] A5755/A6584 
[O n] A7330/A3726 
[Sir] A4068/A6717 
[Oni] A4363/A5007 
[Sin] A6312/A9069 

Average 
We (cm" 3 ) 
T e (K) 



A^cm- 3 ) 

140 ± 30 
150 ± 50 



-50 

100 ± 30 



+ 2000 



510-1450 



9900±200 
9400±900 
10000±500 
10000±200 



(K) 
14800±400 
13100±400 
10000±500 
9800±100 
9500±100 



90 ±40 
10000±300 



.27 



Abundances 

O+/H+(xl0 4 ) 

O 2 +/H+(xl0 4 ) 

O/H(xl0 4 ) 

12+log(0/H) 

S+/H+(xl0 6 ) 

S 2 +/H+(xl0 6 ) 

ICF(S) 

S/H(xl0 6 ) 

12+log(S/H) 

N+/H+(xl0 6 ) 
ICF(N) 
N/H(xl0 6 ) 
12+log(N/H) 



Ar 2 +/H+(xl0 7 ) 
Ar 3 +/H+(xl0 7 ) 
ICF(Ar) 
Ar/H(xl0 7 ) 
12+log(Ar/H) 



0.68±0.09 
1.21±0.17 
1.88±0.36 

+0.08 
-0.09 



150 ± 40 
9800±300 



0.49±0.06 
1.91±0.17 
2.40±0.35 



8.38 



+0.06 
-0.07 



.37 



,50 ±0.02 



0.46±0.03 
2.94±0.15 

1.28 
4.37±0.36 
6.74 4 



,+0.03 
-0.04 



1.97±0.13 

2.78 
5.47±0.37 
6.74±0.03 

7.19±0.36 
0.63±0.11 

1.14 
8.91±1.46 

+0.07 
-0.09 



0.51±0.04 
4.23±0.23 

1.57 
7.47±0.66 
6.87±0.04 

2.10±0.15 

4.92 
10.30±0.73 
7.01±0.03 

6.67±0.51 
0.56±0.07 

1.04 
7.52±1.02 

c oQ+0.06 
O-»«_0.07 



.87 6.99 ±0.10 



7.07 7.21 ±0.08 



6.07 6.26 ±0.10 



a Values taken fromlRussell Donital ll99(Tft . 
6 Values taken from lPeimberj 12003^1 adopting i 2 



6.1 



0.033. 



Empirical technique: line luminosities 

6.1.1 Determining the dominant WR subtypes 

Empirically, the WR population is estimated from av- 
erage line luminosities of WR stars and is thus depen- 
dent on the dominant subtypes assumed to be present 
jSchaerer fc Vaccjll993l . 

For WN stars we use the line width of He II 1640 in 
the UV as the main discriminator of subtype, since it 
can be attributed almost entirely to WN stars; in op- 
tical low resolution data, A4686 is blended with nebu- 
lar emission lines ([Fein] and [Ariv]) or with other WN 
and WC line components (e.g. Nm 4634-41 and Cm /Civ 
A4650/58). Using the generic LMC UV spectra presented in 
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Figure 12. Plots showing th e fits to He n 1640 for cluste rs Al (left panel), A2 (centre panel) and A4 (right panel) using the generic 
LMC UV spectra presented in lCrowther fc Hadfieldl 120061 1 of early (top), mid (middle) and late-type (bottom) WNs. 



ICrowther fc Hadfieldl (2006; hereafter CH06) of early, mid 
and late WNs we reproduce the line morphology of He II 1640 
in Al, A2 and A4 (see Section ETOt . 

In Fig. [H] we show the fits obtained to Hen 1640. The 
apparently good fit from the late-type templates is due 
to the poorer resolution of the template spectra in which 
the apparent He n 1640 line width is due to blending with 
Fe IV 1632. A dominant late-type population can be excluded 
since A1632, if present, would be separated from A1640 at 
the resolution of the STIS UV spectra (3.1 A). Although a 
dominant population of mid-type WNs cannot be ruled out 
entirely, the width of the template line profile is narrower 
than observed in either Al or A2. We therefore conclude 
that the dominant WN population in Al and A2 is early- 
type from the fit to Henl640 3 . For cluster A4 an excellent 
fit to the Hen 1640 line profile is obtained using the mid- 
type templates (see middle right panel of Fig. 1121 and we 
therefore conclude that WN5-6 stars dominate in A4. 

In contrast, for WC stars we use the optical wavelength 
regime to determine the main subtype because C iv 5808 can 
be attributed solely to WC stars, whilst C IV 1550 suffers 
from contamination by O and WN stars. From the presence 
of A5808, and the absence of C m 5696 emission, we conclude 
that WC4-5 stars are the dominant subtype. This is sup- 
ported by comparisons of our measurement of the FWHM 
of the A5808 feature (~ 80 ± 10 A) with those made by 
ICrowther et alJ l|l998t) . We also infer that WC stars are 
solely present in cluster A2 from the strength of C iv 1550 
relative to He II 1640 in the ultraviolet. In summary, we as- 

3 A significant contribution to the He n 1640 emission line flux 
in A2 comes from WC stars and thus the width of A1640 is less 
reliable as a discriminator of subtype in this case. 



sume that early WN (WN2-4) and WC (WC4) stars dom- 
inate the WR populations of Al and A2, with WC4 stars 
absent in Al, while in A4 we assume that mid WN(WN5-6) 
stars dominate. 



6.1.2 Estimating the number of WR stars 

UV: We estimate the WR populations of Al, A2 and A4 
by scaling the generic LMC UV spectra presented in CH06 
to match the intrinsic line morphologies of He II 1640 and 
C iv 1550, as shown in Figs. ll2l and ll3l A comparison of this 
technique versus simple line flux measurements is given in 
CH06. 

For cluster Al we estimate an early WN content of 80 
stars, with WC stars absent. For A2, we estimate the WR 
content by simultaneously fitting Civ 1550 and Hen 1640 
until the strength of He II 1640 is reproduced, although we 
anticipate under-predicting C iv 1550 emission since O stars 
will also contribute to this line (recall Fig. llOtl . We esti- 
mate 10 early WN stars plus ~25 early WC stars in A2. 
Remarkably, in cluster A4 we have been able to detect just 
three mid WN stars despite the large distance of 14.7 Mpc 
to Tol 89, which is testament to the dominance of WR stars 
in integrated (cluster/galaxy) spectra at UV wavelengths. 
To permit direct comparisons with the WR populations in- 
ferred from the optical UVES spectra of knot A, we sum 
the individual WR numbers in Al and A2, as indicated in 
Table we do not include the three WN5-6 stars derived 
for cluster A4 since these will not contribute significantly to 
the optical blue bump which we attribute entirely to WN2-4 
stars. 

Previously, IChandar et alJ l|2004t) derive a value of 
95 ± 68 late-type WN stars based on the their de-reddened 
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Table 5. WR content of Tol89-A derived using the average 
LMC template WR s pectr a and line luminosities (erg s -1 ) of 
lOowther HadfieldlElO^. 





T")i a en n^t c 

Line 


Line Al 
Luminosity 


A 2 


A 1 


A 






Optical 








N(WN2-4) 


A4686 


8.4xl0 35 






100 


N(WC4) 


A5808 


3.3xl0 36 






15 


N(WR) 










145 


N(WC)/N(WN) 










0.5 






UV 








N(WN2-4) 


A1640 


8.4xl0 36 80 


10 




90 


N(WN5-6) 


A1640 


1.7xl0 37 




3 


3 


N(WC4) 


A1550 


2.0xl0 37 


25 




25 


N(WR) 




80 


35 


3 


118 


N(WC)/N(WN) 






2.5 




0.3 



He ii 1640 line flux for Tol 89-1 (which includes our clusters 
Al and A2) and on the average He II 1640 line luminosity of 
a WNL star taken from SV98 (Wo = 1.2 x 10 37 ergss -1 ). 
Prom the line luminosity of A1640 given in their table 3, 
this equates to ~ 130 early WN stars using the CH06 LMC 
line luminosity for WN2-4 stars, a factor of 1.3 times larger 
than our estimate of ~ 100. We have measured the observe d 
Hen line flux from the spectrum of IChandar et all l|2004) . 
kindly made available to us by the author, and have mea- 
sured an observed flux that is in agreement with the sum of 
the observed fluxes in Al and A2. Using the LMC template 
spectra of CH06, we derive an early WN content of ~ 90 in 
agreement with the results obtained for the sum of Al and 
A2. Thus, we attribute this difference in WN numbers (cf 
~ 100 and 130) to the different choice of extinction laws and 
techniques adopted. 

Optical: Applying the same techniques in the optical, we 
estimate the number of WR stars from fits to He II 4686 
(WN+WC) and Civ 5808 (WC) seen in knot A. Recall, no 
broad WR emission features are detected in knot B. Since 
WC stars are known to contribute to the strength of A4686, 
we begin by first estimating the number of WC stars respon- 
sible for the yellow bump (A5808), so that their contribution 
to A4686 can first be taken into account. 

The results of our optical fitting procedure are shown in 
Fig. 1141 From the fit to C IV5808 (top panel, dashed line) we 
estimate a population of 45 early WC stars (cf 25 from the 
UV). Although the early- type empirical fit does not give 
the best match to A4686, a dominant population of mid- 
type WNs can be ruled out from the absence of Niv4058, 
while the UV excludes a late-type dominant population from 
the absence of Feivl632 emission. Taking into account the 
45 early WC stars determined from Civ 5808, we derive a 
population of 100 early-type WN stars, which is in good 
agreement with the results obtained from the UV (cf 90 
WNE for A1+A2), as shown in Tabled 

In summary, our UV and optical WR diagnostics give 
consistent populations of WN stars, plus reasonable agree- 
ment for WC stars. For Tol 89- A (which encompasses Al-4), 
our derived WR content of ~ 95 WNE and ~ 35 WCE stars 




1520 1560 1600 1640 1680 

Wavelength (A) 



Figure 13. Plot showing the fit to C iv 15 50 for cluster A2 using 
the ge neric LMC UV spectra presented in lCrowther fc Hadfieldl 
( 2006) for a mixed population of early-type WN2-4 and WC4 
stars. 



indic ates N(WC)/N(WN) ~ 0.4 (cf ~0.6: ISchaerer et alJ 
1999). 



6.1.3 Estimating the number of O stars 

O- star numbers were derive d following the methods outlined 
in ISchaerer fc Vaccalll99H (hereafter SV98) using the ob- 
served line luminosity of H/3, which gives Qo bs - the total 
number of hydrogen ionising photons. Qo hs is then related 
to the total number of ionising O stars N(O) for a given 
subtype (typically 07V) using Eqn. where the ionising 
contribution from a population of WR stars is taken into 
account: 

_ Qo bs - N(WN)Q1T ~ N(WC)Q^ C 

and where Qp 7V is the Lyman continuum flux of an indi- 
vidual 07V star and r]o(t) represents the IMF-averaged ion- 
ising Lyman continuum luminosity of a ZAMS pop ulation 
norm alised to the output of one equivalent 07V star ijVaccal 
Il994h . N(WN) and N(WC) are the number of WN and WC 
stars respectively. 

From the line luminosities of H/3 given in Table E3 we 
derive the following Qo bs values of ~ 320 and ~ 2000 x 
10 49 s _1 for knots A and B, res pectively. The sum of these 
is ~ 3 times that obtained bv ISchaerer et all lllflflfll ) who 
derived a Qo hs value of 708 x 10 49 s" 1 for the Tol 89 complex. 
This is likely due to the different choice in slit width (l'.'6) 
and PA (39°) (cf Table [TJ. 

Values of r/o (t) are taken from the instantaneous burst 
models of SV98 (see their fig. 21) for the ages given by the 
equivalent width of H/3 obtained in Section II. 2. II We adopt 
r]o(t) values of 0.25 and 0.9 for knots A and B respectively 
(see Table 0. From a calibration of line bla nketed Galac- 
tic O star models, logQ (O7V) = 48.75 s" 1 ^Martins et alJ 
120051) . Recent studies of Magellanic Cloud O stars indicate 
~2-4 kK higher temperatures than their G alactic counter- 
parts jMassev et alJl200.4 Itlean et al.ll2006h . Consequently, 
we adopt IogQ (O7V) = 48.9 s" 1 at LMC metallicity 
jHadfield fc Crowtherl EflOfih . For typical LMC early WN 
and e arly WC star population s we adopt log Qo values of 
49^0 JCrowther fc Sn.i. hi fl99(l and 49.4 jCrowther et al] 
2002), respectively. 



The massive star population in Tol89 15 



45 WC4 




Figure 14. Plots showing the results of our line profile fitting. 
Top: Fit to the C iv 5808 feature using 45 LMC WC4 stars (dashed 
line). 2nd: Plot showing the individual contributions from 100 
early-type WNs (dot-dashed line) and 45 early type WCs (dashed 
line) to the blue bump. The sum of the individual components 
is shown as the solid line. 3rd: As for panel 2 but showing the 
individual contributions from 45 mid-type WNs (dot-dashed line). 
Bottom: As for panel 2 but showing the individual contributions 
from 95 late-type WNs (dot-dashed line). The UVES spectra have 
been binned for clarity, continuum subtracted and velocity and 
extinction corrected. 



Table 6. Nebular derived O star content for Tol 89-A and B based 
on the line luminosity and equivalent width of H/3: L(H/3) and 
W(H/3) respectively. The H/3 line fluxes (F^ and l\) are in units 
of ergss -1 cm -2 . Line luminosities are derived based on D= 14.7 
Mpc and are in units of ergss . The WR numbers for Tol 89-A 
are averages of UV and optical empirical calibrations from Ta- 
ble0 whilst WR stars are not seen in Tol89-B. 



Region 


A 


B 




F(H/3) 


4.7 x 10~ 14 


1.4 x 


io- 13 


I(H/3) 


5.9 x 10" 14 


3.7 x 


io- 13 


L(H/3) 


1.5 x 10 39 


9.5 x 


10 39 


W(H(3){A] 


55 


295 




Age (Myr) 


4.5 


2.5 




/nobs 
qWN 

Qrf c 


3.2 x 10 51 


2.0 x 


10 52 


8.0 x 10 48 


8.0 x 


10 48 


1.0 x 10 49 






2.5 x 10 49 






vo(t) 


0.25 


0.9 




A?WN 


95 






Nwc 


35 






N 


085 


2780 




WR/O 


0.2 







6.2 Synthesis technique: Starburst99 modelling 

In this section we compare the observed optical UVES 
spectrum of Tol 89-A and B to model predictions com- 
puted using the evolut ionary synthesis code Starburst99 v5.0 
l|Leitherer et alJll999F ). The code is an improved version of 
Starburst99, incorporating both a new set of evolutionary 
tracks from the Padova gr oup for old and low-mass stars 
j Vazquez fc Leithererl2005T ). as well as a high-resolution (0.3 
A) optical spectral library (exc luding WR stars) covering the 
full HRD l|Martins et al.ll2005l L 

We have implemented optical WR spectral line synthe- 
sis into Starburst99 by incorporating the high resolution (0.3 
A) University College London (UCL) grids of expanding, 
non-LTE, line-blanketed model a tmospheres for WR stars 
{Smith. Norris fc Crowtherl I2002T ). The UCL models have 
been calculated for the five metallicities represented in Star- 
b urst99 (0.05, 0.2, 0.4, 1 and 2 Z©) using the cmfgen code 
of iHillier fc Millerl l|l99St) . The WR grids assume mass-loss 
scales with metallicity (M — Z ) and replace th e pure helium, 
unblanketed, WR models of ISchmutz et all l|l992l) . Here- 
after, we refer to the computed models as the 'SB99+UCL' 
models. 



For knot A we derive N(0)~ 690 from nebular H/3 
emission, assuming a contribution from 95 WNE and 35 
WCE stars giving N(WR)/N(0) ~ 0.2 (see Table®. For 
region B we derive a content of ~ 2800 O stars from 
Qo bs ~ 2000 x 10 49 s- 1 and 770 (t) determined from the SV98 
models (see their fig. 21). T his is about a f actor two smaller 
than the value obtained bv l.Tohnson et al.1 using radio diag- 
nostics (~ 3800 x lO 49 ^ " 1 , given a dist ance of 14.7 Mpc 
for H = 75 km s- 1 MDC~ 1 : IScha,erer et all) . We expect that 
th e discrepancy b etween our derived value of Qo hB with that 
of I.Tohnson et all is due to the fact that radio observations 
probe deeper into the star-forming region, thus detecting 
massive stars that are otherwise obscured at optical wave- 
lengths. 



6.2.1 Model parameters 

Three models were calculated for an assumed 10 6 Mq in- 
stantaneous burst between and 8 Myr for a 0.1-100 M Q 
Kroupa IMF. We choose to use the Geneva group solar and 
LMC metallicit y stellar evolutionar y tracks with enhanced 
mass-loss rates ( Mevnet et al. Il994|) . For the LMC metallic- 
ity tracks we use both the solar and LMC WR model atmo- 
spheres to test the effects of M — Z scaling (i.e. the default 
scaling with M is turned off). We adopt the Geneva evolu- 
tiona ry tracks over those of the Padova group ijGirardi et alJ 
l200fl and references therein) due to the better treatment of 
assi gning WR atmospheric mod els to the evolutionary tracks 
fsee I Vazquez fc Leithererll2005l) . 



16 F. Sidoli et al. 



6.2.2 Optical 

The massive star content of Tol 89 is obtained by scaling 
the best-fitting SB99+UCL models to match the continuum 
flux of the dereddened spectrum over the wavelength range 
of 3000-7000 A. 

The observed spectrum is first corrected for foreground 
galactic extinction and the necessary amount of internal ex- 
tinction is then applied to the observed data to match the 
slope of the best-fitting model, which is adjusted to match 
the optical continuum flux of the dereddened observed data. 
This scaling, along with the massive star content calculated 
by Starburst99 for our best-fitting (age) model, then gives 
the mass of the burst event, and thus the massive star pop- 
ulation. The predicted N(WR)/N(0) and N(WC)/N(WN) 
number ratios calculated by Starburst99 are, of course, un- 
changed by this mass-scaling. The results of the fitting pro- 
cedure are given in Table 

In Fig. EH we show the best-fitting (age) SB99+UCL 
synthetic spectra to the blue bump in Tol 89-A for each of 
our models. The ages derived are slightly younger than the 
~ 4.5 Myr derived in Section l4~2.il cf 4.0 ± 0.5 (solar) and 
3.5-4.0 Myr (LMC). We can see that the best-fitting model 
to the blue WR bump, and to Hen 4686 in particular, is 
obtained for a 4 Myr solar model; although the solar model 
does not give the best fit to the overall continuum shape. 

While the LMC models do give a better fit to the contin- 
uum, they fail to predict the strength of He II 4686, even with 
M — Z scaling for WR stars switched off. This result is not 
entirely unexpected since the initial m ass for WR formation 
increases with decreasing metallicity l|Mevnet et alJ Il994h 
and therefore we expect fewer WR stars to form for a given 
IMF. In neither of the models is the strength of Nv4620 
matched, although all provide a good fit to Civ 5808. We 
note, however, that even the empirical LMC WR templates 
of CH06 fail to reproduce the observed profile of the blue 
bump (see Section r6.lt . 

We determine the mass of the burst in knot A to be 
in the range (4±1) xlO 5 M© (Solar) and 2-3 xlO 5 M© 
(LMC). Applying the same approach to knot B, we derive 
ages of less than 2.0 and 2.5 Myr, plus masses of 6 and 5 
x10 j Mq, for the Solar and LMC models, respectively. Us- 
ing these masses, and the population predictions made by 
Starburst99, we derive the O and WR star numbers given 
in Table □ Using the LMC models we derive N(O) = 660- 
910, N(WN) = 2-3 and N(WC) = 60-90 for Tol 89-A. The 
N(WR)/N(0) and N(WC)/N(WN) ratios are 0.1 and 17-35, 
respectively. The large N(WC)/N(WN) ratios predicted by 
Starburst99 are discussed in Section0 For knot B we derive 
N(O) = 2030-2120 (Solar) or 1630-1740 (LMC) as shown in 
Table □ 

For Tol 89-A in particular, the O star numbers we de- 
rive from optical continuum fits with Starburst99 are in 
good agreement with the numbers indirectly derived from 
the nebular H/3 line luminosity in Section 16.1.31 cf ~ 690 
with ~ 660-910 for Tol 89-A and ~ 2800 with ~ 1600-1700 
for Tol89-B for L(H/3) versus Starburst99 modelling. In ad- 
dition, the N(WR)/N(0) ratio inferred for Tol 89-A is also 
in excellent agreement, cf 0.2 versus 0.1-0.2, respectively. 
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Figure 15. Plots showing the best fitting (age) SB99 + UCL 
models to the blue and yellow WR bumps of Tol 89-A for solar 
and LMC metallicities. A 3.5 Myr, Z = 0.4Z© metallicity instan- 
taneous burst model for a 0.1-100 Mq Kroupa IMF is shown over 
the spectral range 4400-6000 A. The observed cluster spectrum 
has been velocity corrected and binned for clarity. 



6.2.3 UV 

We apply the same technique to the UV continuum, scaling 
the best-fitting (age) LMC Starburst99 models determined 
in Section 14.1.21 to match the flux levels of the dereddened 
cluster spectra of Al-4 and Bl. We derive the following 
masses for Al and A2 of ~l-2xl0 5 and 2-6 xlO 4 Mq, giv- 
ing N(O) w 310-680 and « 40-190 respectively. For A3 and 
A4 we derive masses of 0.4-1.7xl0 4 and 1-5 xlO 3 Mq, giv- 
ing N(O) « 10-60 and « 5-20 respectively. The sum of the 
masses determined for Al-4 are consistent with the value 
derived from the optical, cf ~l-3xl0 5 and ~ 2-3xl0 5 M Q 
respectively. Similarly, the O star numbers are also in good 
agreement: cf ~ 660-910 and ~ 370-950 from the optical 
and UV modelling respectively. These results are in excellent 
agreement with the ~ 690 O stars derived from the nebular 
H/3 approach. For Bl we derive a mass of ~ 3x 10 4 M Q , giv- 
ing N(O) ~ 80-120. This mass is significantly smaller than 
the ~ 5 x 10 Mq derived from the optical UVES spectrum 
for Tol89-B, which optical STIS imaging reveals is made up 
of at least 5 clusters (see leftmost inset of Fig.~2jl. 
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Table 7. Massive star content derived from SB99+UCL modelling of an instantaneous bursts with a Kroupa IMF (a = 1.3, 2.3) and the 
following mass boundaries: 0.1 (Afi ow ), 0.5 and 100 (M up ) Mq. LMC and Solar metallicity models are presented, in which O and WR 
populations result from optical continuum and emission line fits, respectively. 



Model 






LMC, Z= 0.4Z 







LMC, Z= 0.4Z Q 


Solar 








uv 






Optical 


Optical 


Region 


Al 


A2 


A3 


A4 


Bl 


A B 


A B 


Age (Myr) 


4.5±0.5 


5.0±1.0 


3.0-5.5 


3.0-5.5 


< 3.5 


3.5-4.0 < 2.5 


4.0 ±0.5 <2 


Mass (xlO 5 M ) 


1-2 


0.2-0.6 


0.04-0.17 


0.01-0.05 


0.3 


2-3 5 


4±1 6 


N(O) 


310-680 


40-190 


10-60 


5-20 


80-120 


660-910 1630-1740 


670-1620 2030-2120 


N(WN) 












2-3 


5-20 


N(WC) 












60-90 


90-240 


N(WR)/N(0) 












0.1 


~0.2 


N(WC)/N(WN) 












17-35 


7-31 



7 DISCUSSION 
7.1 Tol 89 as a GHR 

Tol 89 is conspicuous by being the only larg e site of star for- 
mation in NGC 5398. l,Tohnson et 3,1.1^2003) note that Tol 89 
is among the most luminous radio Hn regions yet observed 
and is comparable to NGC 5471 in M 101. 

Tol 89 is composed of three knots of star formation 
which we denote as A, B and C. Optical HST/STIS imag- 
ing resolves the two brightest optical knots of star formation 
(A and B) into multiple clusters. Each is a massive burst of 
star formation in its own right, having produced multiple 
young compact massive clusters (R e ff ^ 3 pc, M >~ 10 3 ~ 5 
Mq) over very short time scales (~ few Myr). The brightest 
cluster in knot A (Al), for example, has a mass of ~1- 
2 x 10 5 Mq , while a total mass of ~ 6 x 10 5 is inferred for 
the ionizing sources at the heart of Tol 89-B. A further two 
young compact (R c ff ~ 2 pc) massive clusters CI and C2 are 
identified in the STIS UV image. The presence of so many 
young compact massive clusters within such an isolated 
GHR makes Tol 89 a rather unusual object. Typical GHRs in 
non-interacting, late-type spiral galaxies tend to host group- 
ings of fewer and less massive (~ 10 3 " 4 Mq) clusters (e.g. 
NGCs5461, 5462 and 5471 in M 101: IChen et a,lJl2005l . and 
NGCs592, 595 and 588 in M 33: IPellerinll2006Jl or multipl e 
OB associations (e.g. NGC 604 in M 33: lHunter et al.lll996>) . 

The fact that Tol 89 is located at the end of the bar in 
NGC 5398 indicates it may have been formed through gas 
inflow in a high pressure environme nt - the conditions un- 
der w hich SSCs are thought to form lElmegreen fc Efremovl 
Il997t) . However. l.Tohnson et a,lJ l)2003F ) suggest that the weak 
bars found in late-type galaxies are not strong enough to 
generate the required gas inflow. It is therefore not clear 
why so many young compact massive clusters have formed 
in Tol 89. 

Larscn fc Richtlerl j 19991. [2000) have shown that mas- 
sive star clusters do form in normal galaxies which show 
no obvious signs of recent interaction - although these typ- 
ically form in isolation. They find that the formation of 
young massive clusters is favoured in environments with ac- 
tive star formation and thus suggest that their formation 
in starbursts or mergers may simply be extreme cases of a 
more general phenomenon. One possible explanation for the 
presence of such a massive star forming region is that the 
parent galaxy NGC 5398 may be undergoing some form of 



interaction. A detailed investigation of NGC 5398 and its 
environment is needed to settle this issue. 

In Table[S]we compare the properties of Tol 89 with the 
three GHRs in M 101 and 30Doradus in the LMC. It can be 
seen that Tol 89 is comparable to, but no more exceptional 
than, the GHRs presented in TableElin terms of its Ha/radio 
luminosities and size. The area normalised star formation 
rate (SFR) of Tol 89 is « 0.1 M yr _1 kpc~ 2 and is com- 
parable to the val ues derived for the cluster complexes in 
M 51 (see table 1 of lBastian et a,ll2005l l and to the definition 
of a starburst galaxy (0.1 Mq yr _1 kpc~ 2 iKennicutt et alJ 
l2005|h The S FR rates in Table Ef are derived using the 
prescription of iKennicuttl l|l998l) : T, S fr (Mq yr _1 ) = 7.9 x 
10- 42 L Ha (ergss" 1 ). 



7.2 Nebular emission lines 

We first consider the origin of the nebular He 11 4686 emis- 
sion in Tol 89-B. This detection is to our knowledge the 
first in a WR galaxy with a metallici ty greater than 0.2 
Zq (cf. iGuseva. Izotov fc Thuanl l20f)ffl . The detection in 
Tol89-A is somewhat ambiguous - due to the poor quality 
of the observations - nevertheless, nebular He 11 4686 is ex- 
pected to originat e from weak-lined 1 metal-poor early-type 
WN or WC stars l|Smith et al.ll2002T) . As for Tol 89-B, the 
absence of WR stars implies that either the most massive 
O stars are responsible for providing the He 11 ionizing pho- 
tons, or that other (non-photoionizing) mechanisms such as 
collisional shocks may be responsible for this emission. 

The intrinsic flux ratio relative to H/3, I (A4686) / I (H (3) , 
is ?s 1 x 10~ 3 . ISchaerer fc Vaccal note that for young 
bursts dominated by O stars (f < 3 Myr), typical val- 
ues for /(He 11 4686)// (H/3) lie between 5 x 10~ 4 and 2 x 
10 -3 (see their fig. 8). Thus, it would seem that our 
/(He 11 4686)//(H/3) ratio is approximately consistent with 



4 The value for Tol 89 is derived from the Ha flux estimated from 
the Ha images obtained as part of the Spitzer Infrared Nearb y 
Galaxies Survey (SINGS) Legacy Project ^Kennicutt et a,lJ2003h . 
The observations were carried out at the Kitt Peak National Ob- 
servatory (KPNO) 2.1 m telescope. We note that this value is 
derived from the observed (i.e. non-extinction correcte d) Ha emis- 
sion fl ux in accordance with the results presented hv lChen et all 
J2005F) in their table 1. 
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Table 8. Comparison of Tol89 with NGCs5461, 5462 and 5 471 in M 101 and 30 D oradus in the LMC. Values for Tol89 are derived from 
Ha images obtained as part of the SINGS Legacy Project iKennicutt et alJ 20031). The observati ons were carried out at the Kitt Peak 
National Observatory (KPNO) 2.1 m telescope. La cul for Tol 89 com es from Johnson et all i2003h . Values for the three GHRs in M 101 
and 30 Dor in the LMC are taken fr om table 1 of lChen et alJ l2005ri . except the star formation rates (SFR) which were calculated using 
the prescription of lKennicutJ il998Tl . The Lhq are not extinction corrected for any the GHRs presented in the table. 



GHR 


Tol 89 


NGC5461 


NGC5462 


NGC5471 


30 Dor 


Angular Size (") 


24 x 18 


40 x 25 


48 x 33 


17 x 17 


1200 x 1200 


Linear Size (pc) 


1700 x 1230 


1400 x 875 


1680 x 1150 


600 x 600 


290 x 290 


^Ha (ergss -1 ) 


1.9 x 10 40 


2.7 x 10 40 


1.3 x 10 4() 


2.2 x 10 4() 


3.9 x 10 39 


SFR Hq (M yr" 1 kpc" 2 ) 


0.07 


0.17 


0.03 


0.47 


0.36 


Lecm 


6.2 x 10 2B 


1.4 x 10 27 


9.1 x 10 26 


7.4 x 10 26 





a young starburst event in the pre-WR phase (in agree- 
ment with the fact the we do not detect any WR emissions). 
However, the SNC02 models - which include a more thor- 
ough treatment of chemistry and line blanketing in their 
WR models - predict I (\4686)/ 1 (H/3) ratios that are a fac- 
tor of fa 10 lower for all ages (~ 1(T 5 and ~ 1(T 4 ), in which 
the lowest metallicity models predict the hardest ionizing 
flux distributions. During the pre-WR phase the SV98 and 
SNC02 models differ due to the use of the costar (SV 98) 
and WM-Basic (SNC02) models fsee lSmith et al]l200l for 
details). 

Consequently, one does not anticipate nebular He II 4686 
in yo ung starburst regions , such as Tol 89-B, to arise from O 
stars. Garnett et discuss Hn regions in which neb- 

ular He II 4686 is observed. In most cases, Wolf-Rayet stars 
provide the ionizing source for these nebulae, except in two 
instances, for which an early O star and a massive X-ray 
binary appear to provide the hard ionization, the former 
being especially puzzling. At present, the source of nebu- 
lar He ii 4686 is unexplained in Tol 89-B, whilst weak-lined 
WN or WC stars could provide the necessary hard ionizing 
photons in Tol 89- A. 

We now consider the origin of the broad velocity com- 
ponents seen in knots A and B. They occur at the same 
velocity which suggests that they have a common origin, al- 
though the knot B component is narrower compared to that 
of knot A (70 vs. 110 km s _1 ). In addition, the broad and 
narrow components appear to have a similar electron den- 
sities. This suggests that the gas responsible for the broad 
component can only be differentiated from the narrow Hn 
region component by its width and that the two distinct 
components probably coexist within the GHR. Finally, we 
find that the maximum velocity of the g measured by 

the FWZI, is 450-600 km s" 1 . 

Underlying broad components to nebular emission lines 
have been reported in a number of studies of GHRs, and 
WR and starburst galaxie s (see iMendez fc Estebanl Il997t 
iHomeier fc Gallagherlll999l . for a review). The properties of 
these components resemble those we have found in Tol 89; 
they have similar ionization conditions to the narrow compo- 
nents, and are spatially ext ended over the star-forming knots 
iMendez fc Estebanl llQQT^I . The origin of the broad com- 
ponent is not well understood; IHomeier fc Gallagher! con- 
sider three possible explanations: (1) it is due to integrat- 
ing over many ionized structures at different velocities; (2) 
it originates from hot, turbulent gas within superbubbles 
created by the winds from clusters; or (3) it is associated 



with some type of break-out phenomenon such as a galactic 
wind. The third option is unlikely because the broad com- 
ponent is not velocity-shifted with respect to the narrow 
Hn region component. The similar densities and ionization 
states of the broad and narrow components argues against 
the second option. We are then left with the first option 
that the broad component is the result of integrating over 
shell structures and filaments at di fferent velocities. Studies 
of the well-res olved GHRs 30 D or flChu fc Kennicutrlll994h 
and NGC 604 l|Yang et al.lll99rJl show that their integrated 
profiles have low intensity broad wings due to fast expand- 
ing shells. We therefore favour this explanation because of 
the similar densities and ionization states of the broad and 
narrow components. 

The narrower width of the broad component in knot 
B compared to knot A may be explained by the younger 
age of knot B (< 3 Myr compared to 4.5 Myr). The cluster 
winds in knot A will be more advanced than those in knot B 
because of the onset of WR winds and supernovae. We would 
thus expect the cluster winds in knot A to have had a larger 
impact on the dynamics and structure of the surrounding 
H II region. This is in accord with the larger reddening we 
deduce for knot B, its morphology showing a higher gas 
concentration (Fig. 4), as well as the fact that we observe 
velocity splitting in the main nebular component for knot A 
but not for knot B. 

7.3 The massive star content of Tol 89 

7.3.1 Clusters Al to A4 

We first consider the difference in the massive star contents 
of clusters Al and A2. While A2 shows a mixed popula- 
tion of WN and WC stars, Al only contains WN stars, al- 
though both have similar ages of 4.5 and 5.0 Myr. A sim- 
ilar difference in the populations of the two subclusters in 
the bright super-star cluster SSC-A in the dwarf starburst 
galaxy NGC 1569 was discussed bv IMaoz. Ho fc Sternberg! 
(2001). Their STIS long slit optical spectroscopy revealed 
the presence of young WR features 5Myr) coexisting 
with an old er red supergiant (RSG) population (Js 4Myr). 
IMaoz et alJ suggest that there is a dichotomy in the pop- 
ulation of the two subclusters NGC 1569-A1 and A2, de- 
spite their similar ages of 5 Myr, in which the WR fea- 
ture originates solely from NGC 1569-A2. They conclude 
that NGC 1569-A1 and A2 must have either widely differ- 
ent IMFs or widely different abundances, or similar, anoma- 
lously high, abundances but slightly different ages. In fact, 
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recent estimates of the age of NGC 1569-A2 have been re- 
vised and it is now believed to be around 12 Myr old 
iAnders et al.ll2004h . thus explaining the difference in the 
stellar populations of NGC 1569-A1 and A2. Of course, in- 
dividual clusters can host mixed (WR and RSG) popula- 
tions if their ages are in the range 4-5Myr. Westerlund 1 
in the Milky Way is an example of a massive cluster whose 
content has be en spatially resolve d into RSG, WN and WC 
populations bv lClark et al.l l)2005h . 

For the case of the clusters Al and A2 in Tol 89 we find 
the difference between Al containing only WN stars, and A2 
both WN and WC stars, can be partially explained in terms 
of age effects due to the rapid evol ution of WR stars on tim e 
scales of typically a few 10 5 years jMaeder fc Mevnetll994) . 
Since the WC phase follows the WN stage for the most mas- 
sive stars, one might suspect that the slightly younger cluster 
(Al) contains only WN stars, while the older cluster (A2) 
contains a mixed population of WN stars and WC stars. 
Unfortunately, one expects the most massive stars to ad- 
vance to the WC stage within ~3Myr, which is inconsis- 
tent with the age inferred from the - Al dominated - H/3 
equivalent width in knot A. Alternatively, the age would 
need to exceed 5Myr for lower initial mass stars to advance 
through to the WN phase, but fail to become WC stars 
prior to core-collapse, which is also somewhat in conflict 
with the observed H/3 equivalent width for knot A. Dif- 
ferences in WR populations could result from variations in 
IMF, with a deficit in very high mass stars for Al, although 
this scenario appears to be rather contrived, given the rapid 
variation in N(WC)/N(W N) with age for young starbursts 
ISchaerer fc Vaccal l|l99SF ). 

A similar situation would appear to arise for clusters 
A3 and A4 in the sense that from their derived ages (see 
Section r4.1.2|l one would also expect WR stars to be present. 
While no obvious WR signatures are detected in cluster A3, 
clear Hen 1640 emission is observed in A4 (see Fig. 0. 
Remarkably, in cluster A4 we have been able to detect just 
three mid WN stars despite the large distance of 14.7 Mpc 
to Tol 89. The detection of such a small number of WR stars 
is partly made possible due to the weak continuum flux of 
this cluster, whose UV output is entirely dominated by these 
three WN stars. In contrast to cluster A4, the observed low 
S/N and stronger continuum in A3 may be masking any WR 
stars present in this cluster, thus we can not exclude the 
presence of a small number WN stars. We derive a mass for 
A3 and A4 of ~ 0.4-1. 7xl0 4 and ~ 1-5 xlO 3 M w from UV 
Starbu rst99 modelling. Accordin g to lCervino fc Mas-Hessd 
jl994 and ICervino et al] j2002|L below ~ 10 4 " 5 M the 
IMF is no longer well sampled, thus affecting the integrated 
properties of the cluster. In such cases stellar population 
synthesis models c an no longer be correctly applied (see e.g. 
I.Tamet et alJ l2004h . It is possible that the massive stellar 
content anomalies of A3 and A4 are due to a stochastic 
sampling of the IMF as a result of the low cluster masses. 

7.3.2 Empirical constraints 

Using template spectra of LMC WN and WC type stars 
JCrowther fc Hadfieldl hoOfih we have been able to derive 
consistent early WN (« 95) and WC (« 35) star numbers 
from both optical and UV diagnostics. In the UV, we sum 
the contributions from clusters Al and A2 in order to make 



direct comparisons with the optical UVES spectra of knot A; 
we do not include the three WN5-6 stars derived for clus- 
ter A4 since these will not contribute significantly to the 
optical blue bump which we attribute entirely to WN2-4 
stars. Based on nebular derived O star populations, we esti- 
mate N(WR)/N(0) ~0.2 for knot A, somewhat larger than 
single star evolutionary models predict at LMC metallicity 
ISchaerer fc Vacc A jl99ST ). 

The UV STIS data are also presented in lChandar et alJ 

i2004h who derive a He 1 1 1640 flux equivalent to ~ 130 early- 
type WN stars (based on the CH06 LMC line luminosity for 
WN2-4 stars) from their extraction of Tol 89-1 (encompass- 
ing our Al and A2). This compares well to the results of this 
work, although we note that the differenc e is most likely due 
to the c hoice of extinction law used by IChandar et al] As 
noted bv lHadfield fc Crowtherl J2006h . the use of a standard 
starburst extinction law is ideally suited to spatially un- 
resolved galaxies. For extragalactic stellar clusters such as 
Tol 89-Al(-A3) and NGC 3125 -A1, an LMC or SMC e xtinc - 
tion law is more appropriate. lHadfield fc Crowtherl l|2006l) 
have shown that the WR conten t of NGC 3125-1 derived 
from He II 1640 bv lChandar et all is strongly overestimated 
due to their choice of extinction law. 

The absolute WR content we derive from the optica l is 
significantly lower than that estimated bv lSchaerer et afl bv 
about a factor of 3; although the WC/WN ratios obtained 
are in good agreement (cf ~ 0.5 and ~ 0.6 respectively). 
This discrepancy is likely due to several factors, including 
different adopted internal extinctions, slit widths (cf l'.'6 to 
our l"4) and p osition angle (PA 39° and 90°). Note that 
ISchaerer e^hj adopt the global extinction value derived for 
Tol 89 bv lTerlevich et"afl l|l99ll) (0.12 mags), whilst we de- 
rive an extinction value of zero mags for Tol 89-A directly. 

7.3.3 Starburst99 model constraints 

In Section l6~2*l we compared three SB99+UCL models to the 
observed WR profiles in knot A; 1) Solar; 2) LMC with 
M - Z for WR stars switched on; and 3) LMC with M-Z 
for WR stars switched off. Fig. El shows that while we are 
able to obtain good fits to the yellow WR bump, the fit to 
the blue bump is rather unsatisfactory - particularly in the 
case of the LMC metallicity models. The best agreement is 
achieved for solar metallicity models, which we acknowledge 
to be unphysical in view of the low metallicity of Tol 89. 

Using the LMC models, the predicted number of WR 
stars is in reasonable agreement with the empirical results, 
approximately a factor of ~ 1-2 times smaller (cf ~ 60-90 
to 130), while the O-star numbers are also similar from the 
stellar continuum (~ 660-910) and Hn region (~ 690) anal- 
yses respectively. The SB99+UCL spectral synthesis models 
predict N(WR)/N(0) ~ 0.1 versus ~ 0.2 from direct stellar 
(WR) and indirect nebular (O) results. Unfortunately, the 
major failure of the spectral synthesis approach relates to 
the distribution of WR stars, which is observed to be pri- 
marily WN stars in contrast to a predicted dominant WC 
population. As a consequence, the fit to the blue bump is 
poor since the WN population is greatly underestimated. 

However, let us recall that non-rotating evolutionary 
mo dels are at present used in Star burst99 synthesis mod- 
els. iMevnet fc Maederl l|2003l . 1200,51) have shown that the 
inclusion of rotational mixing in their evolutionary models 
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increases the WR lifetimes - mainly as a result of the in- 
creased duration of the H-rich phase - and lowers the initial 
mass limit for the formation of WR stars. The combination 
of these two factors is to decrease the N(WC)/N(WN) ratio 
and increase the predicted N(WR)/N(0) ratio, bringing pre- 
dictions closer to our empirical results for Tol 89-A. The po- 
tential use of synthetic WR line bumps in Starburst99 mod- 
els as diagnostics of WR populations is at present severely 
hindered by the lack of rotational mixing in evolutionary 
models. Until then, one should treat detailed N(WR)/N(0) 
and WR subtype distributions predicted by such models 
with caution. 



8 SUMMARY 

We have presented new high spectral resolution VLT/UVES 
spectroscopy and archival HST/STIS imaging and spec- 
troscopy of the giant Hn region Tol 89 in NGC 5398. 

From optical HST/STIS imaging we resolve the two 
brightest optical knots of star formation, Tol 89-A and B, 
into individual, young compact massive clusters Al-4 and 
Bl (Reft < 3 pc, M >~ 10 3 ~ 5 M©). We derive ages for knots 
A and B of ~ 4.5 and < 3 Myr respectively. We determine 
a mass from Starburst99 modelling for the brightest cluster 
in A (Al) of ~ 1-2 xlO 5 M©. A total mass of ~ 6 x 10 5 M 
is inferred for the ionizing sources at the heart of knot B. 
A further two young compact massive clusters (R c ff < 2 pc; 
CI and C2) are identified in the STIS UV image which may 
also fall into the SSC (mass) category. In total, we identify 
at least seven young massive compact clusters in the Tol 89 
star-forming complex. We find that the GHR properties of 
Tol 89 are similar to the three GHRs in M101 and 30 Dor. 
Tol 89, however, contains six clusters of SSC proportions 
whereas the other four comparison GHRs do not. Tol 89 is 
therefore exceptional in terms of its cluster content and its 
isolated location in the late-type ga laxy NGC 5398 . 

In agreement with the results o flSchaerer et all 1999, we 
show that the WR emission is localised to the region with 
maximum stellar continuum, knot A; while STIS UV spec- 
troscopy reveals that the WR stars are confined to clusters 
Al, A2 and A4, with early WC stars located only in A2. 
We have modelled the observed WR line profiles using the 
empirical template spectra of LMC W N and WC stars pre- 
sented in lCrowther fc Hadfieldl l)2006f) . revealing ~ 95 early 
WN stars and ~ 35 WC stars in Tol 89-A. WR populations 
inferred from our empirical technique are consistent between 
optical and UV diagnostics, and so are well constrained. For 
clusters Al, A2 and A4, we obtain N(WC)/N(WN) ~ 0, ~ 3 
and ~ 0, respectively. It is feasible that the slight difference 
in ages of clusters Al and A2 (cf ~ 4.5 and ~ 5 Myr re- 
spectively) is responsible for this difference in their massive 
star population, although differences in IMF cannot be ex- 
cluded. In cluster A4 we have been able to detect three mid 
WN stars despite the large distance to Tol 89, which is testa- 
ment to the dominance of WR stars at UV wavelengths. The 
detection of so few WN stars may be the result of stochastic 
sampling of the IMF. From nebular H/3 emission, we obtain 
N(0)~ 690 and 2800 for knots A and B, from which we infer 
N(WR)/N(0)~ 0.2 for the former region. For knot B, N(O) 
is a factor of two sm aller than recent radio observations by 
I.Tohnson et~aT1 l|2003h . 



We have const ructed complementary Starburst99 
Leitherer et al]ll999l) models in which optical spectral syn- 
thesis of WR stars has bee n imple mented using the UCL WR 
star grids of ISmith et alJ l)2002h . O star populations from 
optical continuum flux distributions are in good agreement 
with nebular results, although the SB99 + UCL models fail 
to reproduce the observed strength of the blue WR bump, 
because too few WN stars are predicted in the evolutionary 
models. The inclusion of evolutionary tr acks with rotational 
mixing jMevnet fc Maederl Eool. 120051) should help to re- 
solve this issue by increasing the lifetime of WNs during the 
H-rich phase. Nevertheless, the total WR populations ob- 
tained in this way agree with the empirical results to within 
a factor of 2-3. 

From an analysis of the optical nebular emission lines, 
we confirm previous determinations that Tol 89 has an LMC- 
type metallicity. Nebular He II 4686 is observed in Tol 89-B 
and perhaps Tol 89-A from UVES and STIS spectroscopy. 
The latter is expected for early-type WR populations at 
low metallicity, whilst He II 4686 emission from starburst 
regions prior to the WR stage is not predicted, unless this 
is formed within a shocked region surrounding young O star 
populations, as may be the case for N44C in the LMC. We 
detect symmetrical broad components in the strongest neb- 
ular lines with velocities up to 450-600 km s" 1 . We find that 
this high velocity gas has similar properties to the Hn gas as 
revealed by the narrow components, and suggest that it is 
the result of integrating over shell structures and filaments 
at different velocities within the GHR. 
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